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Abstract 

We study a particular kind of chaotic dynamics for the planar 3-centre prob- 
lem on small negative energy level sets. We know that chaotic motions exist, 
if we make the assumption that one of the centres is far away from the other 
two (see [7]): this result has been obtained by the use of the Poincare-Melnikov 
theory. Here we change the assumption on the third centre: we do not make any 
hypothesis on its position, and we obtain a perturbation of the 2-centre problem 
by assuming its intensity to be very small. Then, for a dense subset of possible 
positions of the perturbing centre in R 2 , we prove the existence of uniformly hy- 
perbolic invariant sets of periodic and chaotic almost collision orbits by the use 
of a general result of Bolotin and MacKay (see [5]). To apply it, we must 
preliminarily construct chains of collision arcs in a proper way. We succeed in 
doing that by the classical regularisation of the 2-centre problem and the use of 
the periodic orbits of the regularised problem passing through the third centre. 

keywords: Centre problem, Regularisation, Collisions, Chaotic motion. 

1 Introduction 

We consider the motion of a particle in the plane, under the gravitational action of 
three point masses at fixed positions (the planar restricted 3-centre problem). We 
fix a Cartesian reference system Oxy on the plane and choose suitable dimensionless 
coordinates such that the two centres with greater masses occupy the positions C\ = 
(1, 0), C2 = (—1,0). Following the common terminology, we refer to these as primaries. 
We suppose for simplicity that the primaries have equal intensities a\ = a-i = a > 
(symmetric problem), and we always consider all the three centres having positive 
intensities. 

Let C — (xo,yo) eK 2 \ {C\, Ci\ be the position of the third centre and e > be 
its intensity. We assume e to be very small and consider the limit e — > 0: this means 
that £ is a perturbation parameter. In other words, we make the hypothesis that the 
conditions are such that we can consider the problem as a one-parameter perturbation 
of an integrable one: the 2-centre problem. 

Let M be the smooth Riemannian manifold M = M. 2 \ {C\, C2}, with the induced 
Euclidean metric on the tangent bundle TM . The configuration space for the motion 
of a particle in the potential field generated by the centres C±, C%, C, is M\ {C}. Let 
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P = (x, y) denote the position of the particle on M\{C}. Then the system has smooth 
Lagrangian function on T(M \ {C}) given by 



y/(x - x Q ) 2 + (y- y Q ) 2 ' 
with Lq the Lagrangian of the 2-centre problem 



x 2 + y 2 

LO = 7i 



^{x + l) 2 +y 2 y/{x-l) 2 + y 2 



Here we use the Newtonian notation for derivatives with respect to time: x — dx/dt, y — 
dy/dt. Note that Lq is a smooth function on M, while L £ has a Newtonian singular- 
ity at the point C. To simplify notation, let us denote by W(x, y) and eV(x, y) the 
potentials due respectively to the primaries and the third centre C, so that 

W{x,y) = - 
V(x,y) = 



^(x + l) 2 + y 2 ^/( x -l) 2 + y 2 
1 



v /(x-x ) 2 + (y-yo) 2 
The Lagrangian and the Hamiltonian of the problem have the form 

L e = L -eV, H E = Hq + eV , 

with Lq, Hq the Lagrangian and Hamiltonian of the 2-centre problem 

x 2 + y 2 P 2 T + pl 

L = —^-W, H » = F -^pL + W . 

For negative values of the energy H e = E < 0, with E — > 0, we prove the existence 
of chaotic motions passing arbitrarily close to the perturbing centre C. This is made 
by the use of the shadowing result proved by Bolotin and MacKay in [1] . 

For E > non-integrability has been established by Bolotin (see [2]): he proves 
that for the n-centre problem on the plane R 2 , with n > 2, it does not exist an 
analytic integral of motion which is non-constant on the energy shell H~ 1 (E) J with 
E > 0. In [3] the same author extends this result to a wider class of Lagrangian systems 
defined on any 2-dimensional manifold M, with n Newtonian singularities on M: he 
shows non-integrability when n is greater than two times the Euler characteristic 
of the manifold, n > 2x(M), and the energy is over a suitably defined threshold, 
E > E t h- In particular, we have analytic non-integrability for the restricted circular 
many-body problem, in which a particle moves in a rotating plane, under the action of 
the gravitational attraction of n centres fixed on this plane, when n > 2. Nevertheless, 
this generalisation does not add any additional information about the n-centre problem 
on a fixed plane, in which case it reduces to the result given in [2]. 

More recent outcomes exist in literature, valid for positive energy values E > 0: 
see [B], [T3]. In [SJ the authors study the n-centre problem in K 3 , showing that if 
n > 3 and E > 0, then the topological entropy is positive. In [T3], it is proved that no 
analytic independent integral exist for the n-centre problem in the space, if n > 3 and 
the energy is greater than some threshold E > E t h- Moreover, for both the planar and 
the spatial problem, smooth integrals are found on energy levels with E > E t h, where 
E t h = in the planar case. 
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The case E < has been investigated in [7]. The authors study the restricted 
3-centre problem on the plane, when the third centre is very far from the other two 
and consider small negative energies E, in the limit E — > 0. Then they have a two- 
parameter perturbation of the 2-centre problem on the zero-energy level. They succeed 
in applying the Poincare-Melnikov theory, thus proving the existence of chaotic mo- 
tions. 

We too study the case of small negative energies, but our point of view is different, 
in that the position of the third centre does not go to infinity. This allows us to study 
the orbits which undergo close encounters with the perturbing centre C . 

As in [7j the unperturbed system is the 2-centre problem with Lagrangian Lq. After 
the classical regularisation of singularities, obtained by the use of elliptic coordinates 
(£, 4>) £ R x S 1 and a reparametrisation of time, the system separates: we have the 
motion of the standard pendulum for the variable <f> and a two-well potential for £, 
with a hyperbolic equilibrium for £ = 0. We are interested in the motions which go 
very close to the centre C, then we consider trajectories of the unperturbed problem 
that pass through C. The perturbing term of the Hamiltonian is not e-small in the 
vicinity of a collision with C and this results in a lost of control on the stable and 
unstable manifolds, so that a perturbative study of their intersection is not possible. 

Anyway, we will see that for e small enough a hyperbolic invariant set with chaotic 
dynamics exists, for a dense set of possible positions of the third centre in R 2 . More- 
over, the obtained trajectories shadow chains of collision orbits through the third centre 
C. 

We will apply a general result of Bolotin and MacKay (see [1]), valid for Lagrangian 
systems with Newtonian singularities on a d-dimensional Riemannian manifold, with 
d = 2, 3: they show that if the singular part of the Lagrangian is a perturbing term, 
then for each chain formed by arcs of the unperturbed problem which start and end 
at a singular point (collision arcs), there is a unique orbit with the same energy 
which shadows it. Indeed, after local regularisation of the perturbing singularities, 
the corresponding equilibria are hyperbolic and, by the use of the strong A-lemma, 
shadowing arcs are obtained locally in small neighbourhoods of the singularities. It is 
then proved that these arcs can be joined to gain the entire shadowing orbits. In fact, 
the latter correspond to nondegenerate critical points of a formal functional defined on 
small spherical neighbourhoods (circular in the 2-dimensional case) of the singularities. 
This functional can be seen as the energy functional of a Frenkel-Kontorova model. 
There is a correspondence between the equilibrium states of Frenkel-Kontorova models 
and the orbits of symplectic twist maps, and in this correspondence nondegeneracy 
of the equilibrium states of the energy functional corresponds to hyperbolicity, as 
shown in pp. With this approach the existence and hyperbolicity of the shadowing 
orbits is proved and, by proper estimates on the mixed second derivatives of the 
formal functional, it is also derived that the Lyapunov exponents of the corresponding 
Poincare map are of order log e -1 (see [5]). 

The result of [3] has been applied by the two authors to the planar circular restricted 
problem of three bodies, in which a massless particle moves under the gravitational 
attraction of the other two bodies, the primaries, and the latters are supposed on a 
circular orbit about their centre of mass. The second primary is supposed to have 
small mass e. The result is the existence of periodic and chaotic orbits which undergo 
consecutive close encounters with the smaller primary and which approach, as e — > 0, 
arcs of Kepler ellipses around the first primary, starting and ending at collision with 
the second. 

A similar existence result has been obtained also in by a completely differ- 
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ent method. A direct study of the orbits which have close encounters with the small 
primary is carried out and the proof in this case is constructive. Indeed, an approxi- 
mation of the first return map, defined on a region of the phase space whose projection 
is a small circle around the second primary, is explicitly computed. It is shown that 
the first return map is horseshoe like and it allows to conclude about the existence of 
orbits with consecutive infinite close approaches with the small primary. A complete 
numerical study of these orbits is carried out in |12) . 

As previously outlined, in this paper we follow the variational approach of Bolotin 
and MacKay applied to the 3-centre problem: in this manner we don't get an horseshoe 
map, but we still have a symbolic dynamics on the set of quasi-collision orbits found. 

In our case the unperturbed system is simply the 2-centre problem with Lagrangian 
Lq and the perturbing term is given by the potential due to the third centre C. Fixed a 
small negative value of the energy E < 0, by the classical regularisation of the 2-centre 
problem the system separates and we can easily find periodic orbits with energy E 
passing through C . We prove by the implicit function theorem that for almost all the 
possible positions of the perturbing centre C on the plane, these orbits do not collide 
with the primaries C\,C%, if \E\ is chosen sufficiently small. Then they are periodic 
orbits of the not-regularised 2-centre problem and we can use them to define a finite 
set of collision arcs, starting and ending at C, and take as collision chains the infinite 
sequences formed with the chosen arcs. We can verify that the considered orbits of 
the unperturbed problem are nondegenerate (see Section [2] for the definition) and then 
the theorem of Bolotin and MacKay can be applied to obtain a hyperbolic invariant 
set formed by the orbits that shadow the collision chains. Finally, taking as alphabet 
the finite set of collision arcs chosen, a symbolic dynamics is naturally defined on the 
set of shadowing orbits by considering the shift on the collision chains. We can also 
get easily the positivity of entropy. 

In Section [2j we will recall the results of [4] and [5] for the case at hand , starting 
from some basic definitions, and state our main result. In Section [3] we will perform 
the classical regularisation of the 2-centre problem and compute the periods of the 
resulting separated problem. This will allow us to get infinite classes of periodic orbits 
with fixed energy E, passing through the centre C (Section 0]) : this is the core of the 
paper. The delicate point here is to show that the periodic orbits of the regularised 
problem previously found do not meet the primaries for any small enough value of the 
energy. In Section [S] we will verify that these trajectories satisfy the conditions which 
allow to apply the results of [3] , [5] and conclude our proof. 

2 Basic definitions and theorems. 

Consider the unperturbed system with Lagrangian Lq. A solution 7 : [0,T] — > M of 
fixed energy E for this system will be called a collision arc if 7(0) = j(T) — C, and 
7(t) 7^ C for any t £ (0,T). In particular, the latter condition means that there are 
no early collisions. 

Fix an energy value E < 0, such that C belongs to the set D = {(x,y) G 
M\ W(x, y) < E}. Denote by il the set of W 1,2 curves lying in D, starting and 
ending at collision in C. A collision arc 7 of energy E is a critical point of the 
Maupertuis-Jacobi functional Je defined on fl by 
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with 

g E {x,y,x,y) = ^2(E - W(x, y))(x 2 + y 2 ) . 

The arc 7 is said to be nondegenerate if it is a nondegenerate critical point of the 
functional Je on fi. 

This definition of nondegeneracy is the most natural, but it is quite complicated 
to be used for verifications in concrete examples. There are at least four equivalent 
definitions of nondegeneracy, for which we refer the reader to [4]. In particular, a 
sufficient condition for nondegeneracy is the following. Denote the general solution of 
the system (L ) by q(t) = f(qo, fo, t), where qo G M, vq & TM qo are the initial position 
and velocity. Then a collision arc 7 with energy E corresponds to a solution (vq, T) of 
the system of three equations 

f(C,v ,T)=C, H (C,v ) = E. 

If the Jacobian of the system at this solution is nonzero, then 7 is nondegenerate. 
Actually this is the characterisation of nondegeneracy which will be used in Section [5] 
Suppose that the system (Lq) has a finite set of nondegenerate collision arcs through 
the centre C, jk ■ [0,Tfc] — > D, k E K, with the same energy E, where K denote a 
finite set of labels. A sequence (7fc;)i£Z, ki € K , is called a collision chain if it satisfies 
the condition of direction change: 

j kl (T ki ) ± ±7fe, +1 (0) , for any i e Z . 

Collision chains correspond to paths in the graph T with the set of vertices K and the 
set of edges 

T = {(k 1 k')eK 2 \j k (Tk)^±ik'(0)} ■ (1) 
We are going to use the following results: 

Theorem 1 (Bolotin-MacKay, 2000) Given a finite set K of nondegenerate colli- 
sion arcs with the same energy E, there exists £q > such that for all e G (0, £q] and 
any collision chain (7fe i )igZ; ki £ K, there exists a unique (up to a time shift) trajec- 
tory 7:K— > D \ {C} of energy E of system (L £ ), which shadows the chain ('fki)iez 
within orders. More precisely, there exist constants B,B' > 0, independent of e and 
the collision chain, and a sequence of times (ti)ieZ; such that \ti+i — U — T ki \ < Be, 
dist(7(t),7 fci ([0,T fei ])) <Be, for U <t<t l+1 , and dist(~/(t),C) > B'e. 

From Theorem [T] it follows that there is an invariant subset A e on the energy shell 
{H s = E} on which the system (L £ ) is a suspension of a subshift of finite type. The 
subshift is given by the shift on the set of paths of the graph T, that is on the set of 
all the possible collision chains, and the set A e is formed by the orbits which shadow 
them. 

The important fact about the invariant set A e is that it is uniformly hyperbolic. 

Theorem 2 (Bolotin-MacKay, 2006) There exists a cross-section N C {H e = E}, 
such that the corresponding invariant set M s — A e HN of the Poincare map is uniformly 
hyperbolic with Lyapunov exponents of order loge -1 . 

In particular, the set A e is uniformly hyperbolic as a suspension of a hyperbolic 
invariant set with bounded transition times. 

The reader can see [5] for the construction of the cross-section TV and the definition of 
the Poincare map: we omit these technical details here, because it would require some 
machinery, which is not necessary for the sequel. 
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We will show that the assumptions of Theorem Q] are all satisfied for the case at 
hand, for almost all the possible positions of the third centre C in M. More precisely, 
we will prove that, fixed any of these admissible positions of C, for any small enough 
negative energy value, there exist non-degenerate collision arcs, with which we can 
construct collision chains. In particular, our collision arcs are pieces of periodic orbits 
of the 2-centre problem (Lq), which pass through the centre C. We will then derive 
the following 

Theorem 3 Let I C Q + be a finite set of positive rationals. There exists a dense open 
subset Xj C M of possible positions for the third centre C, such that, fixed C G Xj, 
the following is true. There is a small value Eq > 0, depending on C and L, such 
that, fixed an energy value E G (— Eq,Q), we have that there is £q > 0, such that for 
any e G (0, £q) and any sequence (qk)kez, <7fc G /, there exists a trajectory of the planar 
restricted 3-centre problem (L e ) on the energy shell {H e = E}, which avoids collision 
with the third centre C by order e and is within order e a concatenation of pieces of 
periodic orbits for the planar restricted 2-centre problem (Lq), passing through C and 
of classes qk,k G Z (see Subsection \3.3\ for notation). 

The resulting invariant set formed by these orbits is uniformly hyperbolic. 

In particular, as it follows from Theorem [2j fixed a small enough energy value E < 
and e > 0, there is a cross section in the energy shell {H e — E}, such that the as- 
sociated Poincare map has a chaotic invariant set with Lyapunov exponents of order 
log£ _1 and it contains infinitely many periodic orbits, corresponding to periodic col- 
lision chains. The topological entropy of the Poincare map is 0(e)-close to that of 
the topological Markov chain associated to the graph T, defined by (JTJ) . It's easy to 
verify that the finite set of collision arcs that we construct in the proof of Theorem [3] 
determine positive topological entropy (see Appendix |X| . 

Theorem [3] is still true if we substitute in the statement the set Xi with a set X, 
which is dense in M and is independent of the set of rationals /. In this case, we do 
not know if the set X is open or not, the only thing that we can say about it is that 
it is dense (see Remark [6]). 

3 Periodic orbits for the regularised problem. 

In order to apply Theorem [TJ we must find collision arcs with fixed energy E < for 
the unperturbed system (Lq). A natural choice is to look for periodic orbits through 
the third centre C. As a first step, we recall the classical regularisation of singularities 
of Euler. A deep analysis of the 2-centre problem was made by Charlier in the first 
volume of [5]. Following his ideas we get the existence of infinite classes of periodic 
orbits for the separated problem obtained after regularisatiorO]. 

3.1 Regularisation 

We introduce the elliptic coordinates, defined by the map x + iy = cosh(£ + i(f>) from 
the cylinder K x S 1 to K 2 . This transformation has two ramification points at the two 

1 A classification of the periodic orbits of the 2-centre problem based on the variation of the energy 
parameters can be found in [9] and |10| . We will not need this classification: according to it our orbits 
are all of the same kind, because of the constraints that we impose on the parameters. 
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primaries C\ = (0,0), C2 = (0, n). In elliptic coordinates the Lagrangian L e becomes 

Le = ^-^(cosh 2 £ - cos 2 <j>) - W(t 4>) - eV(^ 4>) , 
where the potentials W, V have the form 

unc a\ 2a cosh g 

cosh § — cos z <p 

V(L4>) = 



\J cosh 2 £ — sin 2 <j) + (xq + y 2 ) — 2 (xo cosh £ cos <j) + yo sinh £ sin </>) 

Consider the problem on the energy level set {H e = E}. The Hamiltonian in elliptic 
coordinates is 

•He — E — 5 ^ — H £ > 

cosh £ — cos 2 (j) 



with 



H £ = H - 2a cosh f - [£ - eV(£, 0)] (cosh 2 £ - cos 2 , 



where the symbols p^,p^> denote the conjugate momenta. 

The orbits of the problem with Hamiltonian (H £ ) on the energy level {H £ = E} 
are, up to time parametrisation, orbits of the system with Hamiltonian H £ on the 
energy level {H £ = 0}. The regularised Hamiltonian H £ has no singularities at C\, C^: 
this means that an orbit for H e is an orbit for H e only if it does not pass through the 
primaries. The new time parameter r is given by 



-/ 

Jo 



' ds , (2) 



cosh £(s) — cos 2 (j)(s) 



and this formula allows us to pass from a solution for the system (He), on the zero 
energy level, to a solution for (H £ ) with energy E. 
The Lagrangian corresponding to He is 

C £ = ^ ^ ± ^ ^ + 2a cosh f + [£7 - eV(f, <t>)} (cosh 2 £ - cos 2 0) , 

where the prime sign denote derivation with respect to the new time parameter: £' = 
d£/dr, <jt = d<j>/dT. 

3.2 The separated problem. 

We have to find orbits of the 2-centre problem (L ) with fixed energy E < 0, then 
we put e — and study the regularised system on the energy level {Ho = 0}. The 
Lagrangian is 

(£1)2 1 ( ( k>)2 

Co = 1 ' K ' + 2a cosh^ + £(cosh 2 £ - cos 2 0) . 

The system separates and we have the two one-dimensional problems 

( - 2a cosh ^+ |£|cosh 2 £ = -E x , , 

1 \E\ cos 2 - Si 
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where we have taken into account that E < 0. If we assume \E\ < a, then the potential 
in the variable £ has a local maximum when £ = 0, with negative maximum value 
\E\ — 2a, and two minima ±£ m , defined by cosh(±£ m ) = a/\E\, where the potential 
takes the value — a 2 /\E\. Furthermore, it goes to infinity for £ — > ±oo. Finally, if we 
choose —E\ > \E\ — 2a, i.e. we consider a level over the separatrix, the motion in £ is 
periodic with two inversion points at and £+, defined by the relation:. 



cosh(£-| 



The motion in <f> is simply the motion of the standard pendulum. Then, in the domain 
Ei > 0, we have rotational closed orbits for (f> and there exist action-angle variables. 
The graphs of the potential energies for £ and separately are plotted in Figures Q] 
and [21 when a = 1 and E = —0.5 . 





Figure 1: Potential energy for the variable £ Figure 2: Potential energy for the variable 
when a = 1 and E = —0.5 . when a = 1 and E = —0.5 . 



After these considerations, we make the following assumptions: 
\E\<a, Ei>0, \E\+Ei<2a. 
For future convenience, we scale the energy parameter E\ and define 

A 1 :=E 1 /2a, (3:=\E\/Ei. 
In this manner the conditions (0J) becomes 

A u /3>0, 2/L4i<l, Ai< 

and the regularised system @ takes the form 



(4) 



1 



4a 



f3Ai cosh z £ 
= $A\ cos 2 + Ai 



A! 



(5) 



(0) 



With the assumptions ([5]) on j3, A\, both the one-dimensional motions are periodic, 
with periods Ti,Ta respectively for £, </>, which depend only on the values of f3, A\. If 
the periods T\, have rational ratio T1/T2 € Q, then the corresponding orbit on the 
cylinder KxS 1 is periodic. To find a periodic orbit for the system (Co), corresponding 
to a fixed value of the energy parameter /?, we have to show that there is at least a 
value of A\ for which Ti,Ti have rational ratio. Before doing that we must compute 
the analytical expressions of the periods. 
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Lemma 1 Consider the system with the assumptions Then the one- dimensional 
motions are both periodic with periods T\, T%, for the coordinates £, 4> respectively, given 
by 



2V20- 1 , N 2^/cF T 

Ti = K(k-l) , T 2 = K (n 2 ) , 7) 

where 

K ^ 2-/T=m ' K2 = IT7r K1 ' K2>0 ' 

and K(k) is the elliptic integral of the first type: 



K(k) = cn _1 (°>«) 



dv 



Proof. It is a straightforward computation, starting from the integral expressions of 
t(£),t(0. ' □ 



Remark 1 Observe that lim,,-^- K(k) = +oo and 

i ^ / ^ < 1 
Ki = 1 <^ S 4 l 

In particular, the second condition means that the motion of £ takes place on the 
separatrix energy level. Then, if (3 < 1 the period T\ goes to infinity as A\ — > ^-j-j : 

lim T\ = +oo . 

Correspondingly, for T 2 we have: 

lim T 2 = +00 . 

We are interested to the limit E — > 0, then we can suppose that the parameter Ei 
is greater than \E\. It corresponds to make the hypothesis that (3 < 1. Our definitive 
assumptions are: 

6(0,1), o<a 1<t ^. 

3.3 Periodic orbits. 

The expressions of the periods T\ , T 2 given in Lemma Q] and the Remark Q] give the 
limits 

lim^^i Ti = +00 , lim^i^o T 2 = +00 , 

lim ^-° T i = 7fe K (75) ' lim ^-T^ T * = ^ K ( vS) ■ ^ 

Moreover, from the definition of K(k), we easily conclude that, fixed (3 £ (0, 1), T 2 is 
a strictly decreasing function of A\ G (0, j^a), while T\ is strictly increasing. Indeed, 
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K(k) is a strictly increasing function of k 2 , k 2 docs not depend on A\ and n\ has 
positive derivative with respect to A\ given by 

By these observations, the following is shown. 

Proposition 1 Let f3 £ (0, 1) be fixed. For any positive rational q £ Q + , there exists 
a unique value A x (f3, q) £ (0, j^s) for A x , such that 

qT 1 (f3,A 1 )=T 2 ([3,A 1 ) . (9) 

In particular, the system (0|) has aperiodic solution in correspondence of the value A\, 
with energy E = —2af3A±. □ 



We have seen that for any (3 £ (0,1) and any positive rational q £ Q + , there 
is a periodic orbit for the regularised system (£o) with energy given by the relation 
E = —2af3Ai. We can classify these periodic orbits, identifying each class with the 
rational number q. Then for any fixed value of the parameter /3 £ (0,1), we have 
exactly one value of A\ for each class q £ Q + . Orbits of different classes do not have 
the same energy E; more precisely we have 

Proposition 2 Let (3 £ (0, 1) be fixed. The function Ai(/3, •) : Q + — > (0, j^n), defined 
by the equality (0), is strictly decreasing. Moreover, there exist the limits 

lim A x (I3,q) = — |— , lim A 1 (f3,q)=0 . 

Proof. Hq>q' then qT^fi, A x ) > q'T x (fi, A x ) and we have 

T 3 (P,A 1 (P,q)) = 4r 1 (J3,A 1 (p,q))></T 1 {p,A 1 (P,q))) . 

T2 is strictly decreasing with respect to A x , while T\ is strictly increasing, then 
A 1 (f3,q)<A 1 (P,q'). 

The monotony of A x (f3, •) assures the existence of the limits 

Li = lim Ai{P,q), L 2 = lim Ax(fi,q) . 



Clearly L\,L<z £ 



and from the knowledge of the limits (|5|). we easily obtain 



the desired values for L X ,L2- □ 



It follows that periodic orbits with many "loops" in £ and few in the variable <p 
tend to the separatrix level for <f>; viceversa, orbits with many loops in <p tend to the 
separatrix level for £. In other words, if we increase only for a single variable the 
number of loops before the orbit closes, we will obtain a limit orbit which does not 
close anymore in finite time. For example, take q = m/n, with m, n £ N. Increasing 
the number of loops in £ corresponds to make m larger and consequently A\ smaller. 
The periodic orbit increases the number of "oscillations" in £, while making the same 
number of revolutions in the variable <j>. As m goes to infinity we have that the 
corresponding orbit makes an infinite number of times the same trajectory in £, tending 
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to close, but without being able to reach the limit values cf> = ±7r/2 in a finite time 
interval, in the future and in the past respectively: in particular, it does not complete 
even a single revolution for (f>. Furthermore, the energy E = —2af3Ai tends to zero. 

We conclude that to form easily a finite set of collision arcs with the same energy, 
we should fix q G Q + and look for collision arcs only in the set of periodic orbits of 
the same class q. 

4 Construction of collision arcs. 

In the previous section we have found infinite classes of periodic orbits for the regu- 
larised problem Cq. We would like to use them to construct collision arcs. 

The next step is then to show that, among the periodic orbits of the regularised 
system (Co), there is at least one which passes through the third centre C. Actually 
this is true if the parameter (3 is sufficiently small. We must also verify that the 
obtained orbits are solutions of the not regularised problem (Lq), that is they do not 
pass through the primaries: this is the most delicate point of the proof of Theorem [3] 
Moreover, in Subsection 14.31 we will face the problem of early collisions. 

4.1 Periodic orbits through the third centre. 

Let (£o,0o) £ R x S 1 \ {(0,0), (0,tt)} be fixed elliptic coordinates for the position of 
the third centre C. Then, among the orbits corresponding to the value A\((3, q), surely 
there is one which pass through the centre C, if £o G (£_ (j3, A\), A\)), where £± 
are the inversion points. Note that in Cartesian coordinates this corresponds to say 
that the centre C lies in the region internal to the ellipse defined by the equation £ = £ + . 
By construction, for each j3 G (0, 1), we have A\ G (0, j^s), then lim^o flA\ = and 



1 + A/ 1 - 4/3i? 
lim cosh(£±) = lim - = +oo . 

P^o /j^o 2/3Ai 

Thus we conclude that 

Proposition 3 Fixed a class q G Q + , there exists /3q > 0, such that for any j3 G (0, /3q) 
there is a periodic orbit of system (Cq), associated with the value Ai(/3, q), which passes 
through C , and the coordinate £o is not an inversion point of the corresponding one- 
dimensional motion in £. □ 

Note that the periodic orbits associated with the same value of A\ differ only for 
the sign of the velocities £' , <fi at the centre C. Thus we have exactly two orbits on 
the configuration space RxS 1 : if one has velocity (£q, <j>' ) at C, the other has velocity 
(— £g, <j)' ). The remaining two possibilities give the same orbits, but with the opposite 
direction of motion. Moreover £q ^ 0, because £o is not an inversion point. Then the 
two trajectories corresponding to A\ meet transversely at C on the cylinder RxS 1 '. 

We remark that there is the possibility that the two orbits coincide: it can happen 
when the trajectory has an autointersection at (^o,<Po) before closing. This is a case 
of early collision and it will be treated in Proposition [5] 

The obtained solutions are not yet the collision arcs that we desire. In fact, they 
are orbits for the regularised system (Cq): to be orbits of the 2-centre problem with 
Lagrangian Lq, it's enough they do not pass through the primaries C\,Ci- This is a 
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delicate problem and to face it we will need a general result about the regularity of 
A\ as function of (3, in a neighbourhood of j3 = 0. 

4.2 Avoiding collision with the primaries. 

In this subsection we will obtain the central result of the paper: we will show that 
for almost all the possible positions of the third centre in M x S 1 , the periodic orbits 
through C corresponding to A\ don't collide with the primaries for any sufficiently 
small value of /3. It means that they are solutions of the not-regularised system (Lq) 
and allows us to proceed with the final verifications, in order to apply Theorem [T] 

First of all we would like to recall a useful property which characterises the periodic 
orbits of the regularised 2-centre problem passing through one of the primaries. Note 
that periodic orbits through the primaries always exist because a primary has £ = 0. 

Proposition 4 Let (3 £ (0,1) arbitrarily fixed. Given q € Q + , let to, n € N such 
that q = m/n and (to, n) = 1. Let 7 be a periodic orbit of the system (£q) associated 
with A\{P, q) and suppose that it passes through one of the centres C\, Ci, which have 
elliptic coordinates (0, 0),(0,7r) respectively. 

If n is odd, then the orbit goes through both the primaries in a period, and the 
collisions happens at a time distance of half the period from each other. 

If n is even, then the orbit passes through only one of the primaries and it happens 
two times in a period, at a time distance of half the period. In the configuration space 
R x S 1 the orbit has a transverse self-intersection at the position of the centre. 

Proof. The system (£0) has the form ©. Without loss of generality we can suppose 
that the orbit 7(7-) = (£(r), 4>{t)) passes through one of the centres C±,C2 at time 
t = 0. The orbit 7 collides with a primary at time r 7^ if and only if £(r) = 
and 0(t) 6 {0, n}. Then we must have r = fc-y- = j-y-, with k,j € Z\ {0}. Then 
h = 21 = q and this implies that there is i € Z such that k = im and j — in, because 
(to, n) = 1. It follows that j(t) is a primary if and only if r = im^- = in^- = i^, 
where T = mT\ = nT^ is the period of the orbit. This concludes the proof. □ 

If we assume q = 1, then, using the above Proposition, we can exclude at once 
the collision with the primaries for some simple cases. To simplify notation we place 
i 1 (/?) = i 1 (/3,l). 

Proposition 5 Let 7 be a periodic orbit through the third centre C , corresponding 
to the value Ai(fi), with f3 < /3q, as in Proposition \Q Let (^Oj'/'o) be fixed elliptic 
coordinates for the centre C. If <fio = k^,k £ Z, or £0 = 0, then the orbit 7 cannot 
pass through the primaries. 

Remark 2 Note that the positions with 4>q = k^,k £ Z, and the ones with £0 = 
correspond in Cartesian coordinates (x, y) to points on the coordinate axes. In partic- 
ular, if £0 = then the three centres are collinear on the x-axis, and C lies between 
the primaries. If 4>q = kit then the three centres are still aligned on the x-axis, but C 
is external. Finally, if 4>o = (2k + then the centre C lies on the y-axis and the 
configuration of the centres is symmetric with respect to this axis. 

In particular, Proposition [3| says that, for q = 1, the periodic orbits through the 
primaries intersect the coordinate axes only at the primaries and for £ = £±. 
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Proof. We are in the case q = 1, then from Proposition 2] we know that the orbit 7 
passes through a primary if and only if for any time r such that £(t) = 0, we have 
4>(t) = kir, k G Z, and viceversa. The centre C does not coincide with a primary, then 
it cannot happen that </>o = kir, k G Z, and £0 = at the same time and in these cases 
the statement is obvious. Now suppose 4>o = (2fc + The time intervals from C to 
a position with <p = in, i G Z, are (2j + 1)^-, J € Z, where T = Tj = T2 is the period of 
7. Look at the variable £: the only positions which have a time distance of (2j + l)-j 
from £ = are the inversion points £± , but £0 7^ £± • D 



Remark 3 iVo£e that the proof of Proposition^ cannot be generalised to arbitrary fixed 
values of q G Q + . For example, if q — 1/2, then, if a periodic orbit passes through C\, 
it certainly passes through a point with elliptic coordinates (£,7r/2), with £ G (0, £+): 
in fact the time passed from the last passage through C\ to this point is T 2 /A = Ti/8. 
Nevertheless, it's easy to see that even for this case the passage through the primaries 
is excluded when £0 — or (fio — kit: indeed, the time between two positions with £ = 
is Ti/2 = T2, and the points with <po — kir along an orbit which passes through a 
primary must have £0 = 0,£±. 

Our next step is to study the regularity of the function Ai(/3,q) with respect to 
the real parameter j3, to be used for proving the central result of the subsection: in 
particular, we are interested in the behaviour near = 0. 

Lemma 2 Let q £ Q + be fixed. Then A\ is a smooth function of (3 € (0, 1) and it can 
be smoothly extended to ft = 0. In particular, there exists the limit 

ii(0,g) :=Hmii(^?)e(0,l), 

the periods T± and Ti are smooth for f3 = and 

qT 1 (0,A 1 (0,q))=T 2 (0,A 1 (0,q)) . 
Proof. Denote by F the function 

FOMx) := [qT 1 -T 2 ](0, A t ) , 

defined on the domain V = Uf3,Ai) €t 2 | /? € (0, 1) , < Ai < ^ j. 

We observe that «i, k 2 are C°° functions of (f3, Ax) £ T), and K\, k 2 € (0, 1); then, 
by derivation under the integral sign, we conclude that F is C°° on the same domain. 
Furthermore J^j- > 0, then from the implicit function theorem we can assert the 

regularity of A\ on (0, 1). 

We want to extend the definition of the function Ay to /3 = 0. First of all we 
prove that the function F can be smoothly extended to (3 = 0. We study T\ and T 2 
separately. We start from T x . If /3 € [0, 1) and A\ < jj^, then (1 - AfiA\) > 0. It 
follows that k\ is well defined and C°° for f3 = and that fcf G (|, 1). Then Ti is C°° 
for /3 = 0. Now see T 2 . For /3 G [0, 1), we have n\ G [0, |) and it is a C°° function of 
/3. To assert the regularity of T 2 , it was enough to show k 2 < 1, then we surely have 
T 2 defined and smooth for j3 = 0. 

The next step is to verify that for j3 = there is a unique value Ai(0,q) G (0, 1), 
such that F(0,ii(0,g)) = and that for this value J^(0, ii(0, q)) > 0. From the 
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definition, we easily see that lim K ^oK(K) = ?, while lim,^!- K(k) = +00. Moreover 



q T 1 (0,A 1 ) = ^K(J^±±), T 2 {Q,A{) 



2a V 2 " v ' ^Ai 

Ti(0, A\) is a strictly increasing function of A\, while T2(0, Ai) is strictly decreasing, 
and 

lim Al ^ F(0, Ax) = -00 , 
lim Al ^ 1 -F(0,A 1 ) = ^= 



We conclude that Ai(0, (7) is uniquely determined from the equality F(0, A\{Q, q)) = 0, 
and Ai(0 : q) G (0, 1). Moreover ^-(0, ii(0, q)) > 0, then the regularity of the func- 
tion Ai in j3 — follows from the implicit function theorem. □ 

Now we are ready to state and show the main result. 

Theorem 4 Let q G Q + a fixed positive rational number. There is a dense open 
subset X' q C K x S , such that fixed (Ccb'Ao) G X' q , there is /3q > 0, such that for 
each f3 G (0,/3o), the periodic orbits through (£,o,4>o), associated with Ai{j3,q), do 
not pass through the primaries C\,Ci- In particular, after scaling time, they are 
orbits of the not-regularised 2-centre problem with Lagrangian Lq, and they have energy 
E = -2a^ii. 

Proof. Let (£0, 4>o) 6 R x 5 1 be the position of the centre C in elliptic coordinates and 
let 7(7-) = (£(t), 0(t)) be a periodic orbit associated with Ai(/3, q), which pass through 
(£o,<^o) with velocity (£ , 4>'q)- Without loss of generality we can assume 4>' > (see 
Subsection 14. 1[) . Suppose that the orbit 7 passes through a primary C\ or C 2 : at that 
instant we have £ = and <fi G {0,7r}. Let q = — , with m,n £ Z positive integers 
and (m,n) — 1. Let At be the shortest time interval to go from a primary to the 
centre C along the orbit 7. Thanks to Proposition 01 we must have At < -j, where 
T = mTi = nTi is the period of 7. 

We have two possibilities, corresponding to start from C\ or from C 2 - The orbit 
solves the separated system ([5]). Then, in the first case 

1 f 00 1 
At = —= =dcf> + iT 2 , 

Wa Jo cog2 ^ + 

while in the second case 

1 Z"^ 1 T 

Ar=— =/ =d4>±-l+i'T 2 , 

WaJo ^jpXx cos 2 + ii ^ 
with i, i' G Z, z, i' > 0, and in the second case with sign ± respectively when <f>Q < 7r, 

The condition At < — = ^f^ implies that i < f and z' < ^y^-- Look now at the 
variable £. In both cases we must have 

1 /" So 1 Ti 

At = ±—^ - m + j— 

V« Jo / cosh £ _ ^ cosh 2 ^ _ ^ 2 
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with j e Z, j > 0, and the condition At < i" = implies that j < m. 
We denote by P(<j) ,(3) and Q{£,o,f3) the functions 



Pi-.,,, n = — / = 1 =d<j>, 

' PAi cos 2 + Ai 



2Va 



Q(Co,/3) = ^ / 1 ^ 



2^ 



cosh£ - (3 A x cosh 2 £ - A x 



From Lemma [2] and the fact that £o is not an inversion point, we deduce that these 
functions are smooth in a neighbourhood of (3 = 0. The condition to pass through a 
primary for the first case is 

At = P(0 O , P) + = ±Q(£ , P) + j y > 

while for the second case it is 

At = P(&, /?) ± y + *'T 2 - ±Q(&, P) + j y . 

We have a finite set of possible values for the integers and the equality qT\ = 

^Ti = T 2 holds, then we can summarise all the conditions with the following one: if 
the orbit 7 passes through a primary, then G + (£o, <fro, ft) € S or G~ (£0, </>0 7 P) G S 1 , 
where S C Q is a finite subset of rationals, depending only on the fixed parameter 
q € Q + , and the functions G + , G~ arc defined by 

The functions G ± are smooth in all the variables: in particular, if G ± (^o, ^0; 0) ^ 5, 
there is /3 > 0, such that G^fa, </>o, P) ^ S 1 for any (3 £ (0,/?o)- Then, to exclude 
the collision with the primaries for small values of /3, it is enough that the centre G 
belongs to the set 

X' q = {(^Ao)eRxS 1 \ G ± (Co,^o,0)^5} . 

The set S is finite, then it is clear that X' q is open and dense in E x S 1 . Indeed, we 
easily see from the definitions of the functions that their partial derivatives with 
respect to <f>o arc different from zero; then X' q is a finite intersection of open dense 
subset of R x S . Moreover, the complement of X' q has zero Lebesgue measure. □ 



Remark 4 When q = 1, it follows easily from Proposition^ that at least one of the 
two periodic orbits through C corresponding to the same Ai(/3) does not collide with 
the primaries. Indeed, without loss of generality we can assume cj>' > (see the end 
of Subsection \4-l\ l- Suppose <pQ £ (0, ^): then the time to pass from the centre C\ to 
the centre C must be less than -j , where T — T\ = T 2 is the period. If £0 > the only 
possibility to collide with the primaries is that £' Q > 0, while if £0 < we must have 
£0 < 0. For the other possible intervals of values of 4>o a- similar reasoning works. 

For our purposes, the fact of having at least one orbit through C that does not 
collide with the primaries is not enough: indeed, after we have obtained the collision 
arcs, we want to construct collision chains with them, and to do it we need at least 
two orbits of the system (Lq) with the same energy E, which pass through C with 
not-parallel tangent fields. The latter condition will be investigated in the next section. 
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Remark 5 Theorem [7] would be improved if we showed one of the following: 

i) the partial derivatives a ^ {£,o,<Po,0) are zero only for isolated values o/(£o;<^o)/ 

ii) the partial derivatives (£oi0o>O) do not vanish for every £o G M., <fio 6 5 1 , 
tmift (^Oj ^0 ) a primary. 

For case (i) we would have that the collision with the primaries is possible only for 
isolated positions (£o>0o) G R x S 1 , while for case (ii) the collision would be excluded 
for any position of the centre C . 

At present, we don't have any of these improvements. 

Corollary 1 Let I C Q + be a finite set of positive rationals. There is a dense open 
subset X'j C K x S 1 , such that for each (£o,</>o) G X'j, there is 0o > 0, such that for 
any (3 G (0,/?o), the periodic orbits through (£o> </>()) corresponding to Ai(j3,q), with 
q G I, do not pass through the primaries. In particular, they are periodic orbits for the 
system (Lq). 

Proof. It is an immediate consequence of the construction of the subset X' q cKx S 1 
in the proof of Theorem 2) we can apply this theorem for any q G /, then take the 
intersection of the sets X' q thus obtained, and the resulting set maintains the same 
properties of the sets X' . □ 



Remark 6 Denote by tp : M x S 1 -> R 2 the map to pass from elliptic to Cartesian 
coordinates, V>(£> <P) — cosh(£ + icj>) = x + iy. The map tp is open and surjective then, 
fixed a finite subset I C Q + , the image Xj = ip(Xj) is an open dense subset o/R 2 . 
Then we can say that the thesis of Corollary [7] holds for any position of the centre C 
in an open dense subset Xi C M. 

Consider now the set X' = ri qe q+X' q : by the classical Baire's Category Theorem 
the set X' is dense mix S . Then, the set X — ip(X') is dense in M 2 and for 
any finite subset I C Q_ + , it is contained in Xi . Then we can say that the thesis of 
Corollary [7] holds for any position C G X C M. In this manner the only information 
that we have lost is that we don't know if the set X is open, while we have gained the 
independence of the dense set X of positions of C from the set of rationals I. Note 
that the choice of a sufficiently small f3 cannot be independent of I , instead. 

4.3 Early collisions 

We defined a collision arc in Section [5] to be a critical point of the Maupertuis- Jacobi 
functional which starts and ends at collision with the centre C and does not meet this 
centre at intermediate times. If we take one of the periodic orbits through C found 
in the preceding section, and we pass to Cartesian coordinates, then we cannot be 
sure that the orbit does not pass newly through C before a period has passed. This 
can happen in two ways: when the orbit, considered in elliptic coordinates, meets the 
point (£o, 4>o) and when it meets (— £o> —4>o)- 

When a periodic orbit starting from the centre C passes newly through C in a time 
shorter than its period, then we talk of early collision. To understand when an early 
collision occurs we need to study the behaviour of periodic orbits a little deeper. 

An easy example is given by the case q = 1: in this case any periodic orbit crosses 
the y-axis, is symmetric with respect to this axis and, if it does not collide with the 
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primaries, it has a transverse autointersection at the point of crossing. The symmetry 
comes from the fact that in a time interval of half the orbit's period we have the 
passage from a point (£, <p) to a point (—£,</> + 7r), that is from a point (x,y) = 
(cosh(£) cos(0), sinh(£) sin(0)), to its symmetric (— x, y). The intersections with the 
y-axis occur when (j) — ±7r/2 and from one intersection to the next there is a time 
interval of half the period. In elliptic coordinates the orbit passes through a point 
(£, 7r/2) and after half a period it arrives at (— £, —ir/2), but these two points coincide 
in Cartesian coordinates and are on the y-axis. Then each orbit autointersects at a 
point on the y-axis. The transformation of the velocities when we pass from elliptic 
to Cartesian coordinates is given by: 

TTit- ±\ f €' \ TTit- x\ ( sinh£cos</> - cosh £ sin </> \ , 
v = Ufa*) • ), Ufa*) = ^ coshesin0 sinHcos ^ ) . (10) 

where £/(£, 4>) is an invertible matrix, except when (£, </>) € {(0, 0), (0, %)}. We observe 
that {/(— £, —0) = —U(£, (j)). If at the point of intersection with the y-axis £ = ±£+, 
then £' = and the velocity has zero y-component: then the two crossings are not 
transverse. But this is the case in which the orbit collides with the primaries, reversing 
its direction at the collisions. If £ ^ ±£+ instead, the velocities at the point of 
intersection with the y-axis have both the components different from zero, then they 
are transverse, because by symmetry they differ only in the sign of their x-componcnt. 

An illustration of this situation is given in Figure [3l where we have drawn 18 
periodic orbits corresponding to the value of = 1/7, included the two orbits that 
collide with the primaries: these latter orbits are in green, while the primaries are 
marked with red asterisks. We have also drawn the periodic orbit £ = £+: we cannot 
use it to construct collision arcs, because it has different energy, but this orbit is 
interesting in itself because it encloses all the trajectories with the same values of the 
parameters (and then the same energy), and it is tangent to all. The behaviour of a 
couple of periodic orbits through the same point can be seen in Figure 21 where the 
starting point is marked with a black asterisk. 




Figure 3: Periodic orbits in Cartesian coor- 
dinates when a = 1, q = 1, P = 1/7 and the 
orbit through the primaries. 



Figure 4: Periodic orbits through (£o,0o) = 
(|£+,0) in Cartesian coordinates, when a = 
1,5 = 1,0 = 1/7. 



Then, if q = 1, we can say that there is an early collision at C if and only if C lies 
on the y-axis: in this case there is a unique, up to reversing the direction of motion, 
periodic orbit through C. 

In particular, we see that the period does not change in the passage from ellip- 
tic to Cartesian coordinates. Actually, this holds for any value of the parameter 
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q € Q + . In fact, if we take a periodic orbit that starts from a point (£o,4>o), which 
is not a primary, we see that the matrix U(^o,<po) given in (1101) is invertible and 
U(— £o> — <^o) = — U(£o,(t>o)- Then the only possibility to have a shorter period when 
we pass to Cartesian coordinates is that the orbit passes through (— £o,~0o), with 
velocity (— £ ,— 4>'o)- This is clearly impossible because the sign of <j>' never changes 
along the orbit. 

Anyway, the situation complicates if q ^ 1 and we cannot get a global view easily: 
as an example, in Figure [5] we have drawn the orbits through the point (£o:0o) = 
(|£+,0), when a = 1, [3 = 1/7 and q — 2. To better see the autointersections we have 
put an enlargement in Figure [6l 
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Fi gure 5: Periodic orbits through (£o,<^o) — 
(|£+,0) in Cartesian coordinates, when a — 
l,q = 2,p = 1/7. 




■2 2 



Figure 6: An enlargement of Figure [5] where 
the autointersections of the two orbits are vis- 
ible. 



We note the following remarkable fact: 

Proposition 6 Let (£o,(fro) G RxS 1 be fixed and consider aperiodic orbit (£(t),4>(t)) 
through this point, corresponding to some fixed values of j3,q. If (£(t),0(t)) passes 
newly through (£o, 4>o) before closing, then it is the unique periodic orbit through (£o, 4>o) 
corresponding to the same fixed values of f3, q, up to reverse the direction of motion. 
Moreover, if £o * s not an inversion point, the orbit has a transverse autointersection 
at (Co, <h)- 

Proof. If the initial velocity in (£o, </>o) is (£cu 4>'o)j then the only other possible velocity 
through the same point along the same orbit is (— £q, (f>' ). If the orbit arrives at (£cb </>o) 
before that a period is passed, then it must happen with velocity (— £o)$))> an< ^ ^ n ^ s 
implies that the two possible periodic orbits through (£o><^o) coincide. □ 

At this stage, we know completely early collisions only for the case q = 1. In 
general, we can only say that early collisions cannot be excluded: we don't have any 
knowledge about the conditions that determine them. Anyway, it is not a problem for 
the proof of Theorem [3] When an early collision occurs, we will take as collision arc 
the partial arc of the periodic orbit through C, which starts from C and ends at the 
first next passage through C. The final time T of the collision arc in this case will not 
be the period, but it will be the time of the first return to C . 

Finally, we want to stress that a collision arc must end at C with one of the 
four possible velocities determined by the choice of the parameters, which in elliptic 
coordinates are {±(£ > 4 > 'o)i i( — Coi $))}> wriere (£o> $>) ^ s * ne initial velocity of the arc. 
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5 Nondegeneracy and direction change. 



In this section we will verify that the collision arcs obtained from Corollary [JJ satisfy 
the nondegeneracy condition and moreover they meet transversely at the centre C: 
this will allow us to apply Theorem [JJ and derive our final result, Theorem [3] 

5.1 Nondegeneracy 

Let q G Q + be fixed and m,n € Z be positive coprime integers such that q = m/n. 
Denote by (£o ; <j>o) € Kx S 1 the position of the centre C in elliptic coordinates. Suppose 
that (£ 0) </>o) € X', where the set X' C M x S 1 is the one given by TheoremSJ Take /3 
small enough, so that the periodic orbits of system (Co) corresponding to Ai(/3, q) and 
passing through C do not pass through the primaries. Let j(t) be one of the resulting 
collision arcs for the system (Lq), which starts from the centre C with velocity vq and 
whose energy is E = —2a(3A 1 ((3, q). Then 7(f) solves the system 

f(C,v ,T) = C, H Q (C,v Q )=E. (11) 

As reminded in Section^ to show the nondegeneracy of 7 it is sufficient to verify that 
the Jacobian of system (fTTj) is nonzero. Actually we will verify a slight variant of this 
condition. 

It is convenient to consider as variables the parameters ((3,Ai) instead of the co- 
ordinates of the initial velocity wq. This procedure is right only if there is a local 
diffcomorphism which allows to pass from vo to (/3,A\). The transformation of the 
velocities in the passage from elliptic to Cartesian coordinates is given by the invertiblc 
matrix U(£o, 4>o) defined in (fTOj) . We know that the orbit 7 corresponds to a solution of 
the separated system ([B]), and in elliptic coordinates we have (£o,4>o) = t?|(0)(£oj $>)> 
with ^(0) = (cosh 2 (£ ) - cos 2 (0o)) _1 . Then v is a C°° function of (££,,<#,). From 
the system ((SJ) we see that (£' , (j>' ) is a C°° function of (f3, A%) and then vq is. By the 
same reasoning we see that (/?, A\) is locally a C°° function of the velocity (£q, (p' ) and 
then of Vq- 

For all nearby trajectories from the same initial point C, we evaluate the m-th 
positive instant of time at which they meet the ellipse £ = £oj with the velocity £ equal 
to the initial velocity (0 ^ 0; then we consider the time distance from the n-th passage 
through 4>o, which is given by nF(j3, A\), where 

F(/3,A 1 ) = [qT 1 -T 2 ](p,A 1 ) . 

This is equivalent to consider the value of the angle coordinate <f> at the instant at 
which we have the m-th oriented crossing of the line £ = £o- 

Note that in this manner the time variable T is fixed as function of (/3,A\), T = 
mTi(f3, Ai), so that we have reduced the order of system (fTTj) of a unit. The orbit 7 
satisfies 

F(/9,Ai)=0, ~2apA 1 = E, 

and it is nondegenerate if the Jacobian determinant of this system at the solution 
(/3 , Ai(/3 , q)) , corresponding to the orbit 7, is different from zero. 
The Jacobian matrix is 

\ -2aA 1 -2a/3 J 
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As observed in Subsection 13.31 we have dF/dAi > 0. Furthermore, we have that 
Ai(j3,q) € (0, jq-g)- These properties remain true for (3 = 0, as proved in Lemma [21 
and moreover T\,T^ are C°° functions of (/3,Ai). It follows that when /? = and 
Ai = A\ (0, q) , the determinant is well defined and different from zero, and by regularity 
the same is true also for values /3, Ai(/3, q), with (3 > sufficiently small. 
We conclude that the collision arc 7 is nondegenerate for /3 small enough. 

5.2 Direction change 

Let q € Q + be fixed and m, n € Z be positive coprime integers such that (7 = m/n. 
In order to construct collision chains to which apply Theorem [TJ we should have at 
least two collision arcs which start and arrive at the centre C with transverse tangent 
fields. 

From the results of Section HJ if in elliptic coordinates our centre (£o:<^o) belongs 
to the dense open subset X' q C K x S 1 , then for small enough values of f3, after 

reparametrisation of time and change of coordinates, the orbits associated to Ai(/3,q) 
do not collide with the primaries and they are periodic orbits for the not-regularised 
two-centre problem (Lo). As observed in Subsection 14.11 there are exactly two or- 
bits corresponding to Ai(/3, q), which pass through the centre (^oi^o) with velocities 
(£q, </>q), (— £q, 4>' ) respectively. Changing sign to cj>' , we obtain simply the same orbits, 
with reversed direction of motion. Moreover, the two transverse orbits coincide in the 
case of an autointersection at (£0 <j>o), as we have seen in Proposition |BJ 

Now consider the two orbits associated with Ai(j3,q): suppose they pass through 
the centre (£o,<A)) at time r = and that their velocities are (£q,0q), (— £o>0o)- 
They are obviously transverse at the point (£o,<po) m the cylinder Kx S 1 , because 
we have chosen j3 small enough to have £q ^ (see Proposition [3]). We must verify 
that this transversality is conserved after time reparametrisation and changing from 
elliptic to Cartesian coordinates. The reparametrisation of time is given by formula 
it maintains the directions, because the centre (£o>0o) is not a primary. As seen 
in Subsection |431 the passage to Cartesian coordinates is given by an invertible matrix 
U (see (JTOJ) ) , then the transversality is conserved. 

Look now at the other possible elliptic coordinates for C: they are (— £0,— <po), 
then the possible velocities at this point are the same as the ones at (£0, </>o)- It follows 
that there are no more orbits through C that we can consider, for the same values of 
/?, q. We conclude that there are two transverse directions for the collision arcs starting 
from C . In correspondence of the same values of the parameters /3, q we obtain four 
collision arcs, divided in pairs of arcs with transverse initial velocities. 

Finally, we observe that, in correspondence of some fixed values of j3, q, if an early 
collision occurs at C, then there is a unique (up to reverse the direction of motion) 
periodic orbit through C with transverse autointersection at C. Indeed, since the 
passage through the primaries is excluded, the matrix U is invertible and we cannot 
have inversion points along the orbit, because cj>' =/= 0. 

We summarise the results of this and the preceding subsection in the following 

Proposition 7 Let I C Q_ + be a finite set of positive rationals. Suppose that the cen- 
tre C has elliptic coordinates {£o,4>o) € X'j, where X'j ClxS 1 is the dense subset 
given by Corollary\]\ Then, there is /3q > such that for each f) € (0,/3o) and for each 
q G / ; there are exactly four collision arcs at C for the system (Lo), associated with 
the value Ai(J3,q), which are nondegenerate. Moreover, the four collision arcs divide 
into two pairs, according to their initial velocities at C: any couple is formed by two 
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arcs with opposite initial velocities at C , and each arc in one pair has transverse initial 
velocity to each arc in the other. □ 



5.3 Proof of Theorem [3] 

We're going to conclude the proof of Theorem [3J 

Let / C Q + be a finite set of positive rationals and Xj = ip(X'j) C M as in 
Remark [51 The energy of the collision arcs at C is given as function of the parameters 

q), q € /, by the relation 

E = -2a/3Ai(/3,q) . 

Thanks to the regularity of the function Ai((3,q) for f3 £ [0, 1) and to the fact that 
Ai(0,q) £ (0,1) (Lemma [2]), to require (3 small enough is equivalent to ask for the 
absolute value of the energy E to be small enough. In particular, if /? < /3o> with /?o 
sufficiently small, then the energy E(f3, q) is a strictly decreasing function of (3: 

£(•,?) :(0,A>) ->(-■$),()), 

with E(0,q) = 0, E(f3o,q) — —Eq < 0. It follows that all the previous results can be 
stated using the energy parameter E instead of j3. 

If we choose the energy E < sufficiently close to zero, then for each q £ I, there 
are four nondegenerate collision arcs of energy E through the centre C: they are pieces 
of periodic trajectories on the configuration space M = M 2 \ {Ci, C2}. 

For f3 fixed, the energy increase with the class q (see Proposition [5]) . This means 
that we cannot state a general result valid for a fixed energy E and any q £ Q + . On 
the other hand, we don't need such a result, because to apply Theorem Q] we only 
want a finite number of collision arcs. What is certainly true is that for any finite set 
of classes I c Q + of cardinality i, we can choose a small enough energy value E < 
to form a set of Ai collision arcs of energy E, by taking for each q £ I the four arcs 
obtained by our procedure. 

By the monotony of the function Ax{q,fi) with respect to q and its relation with 
the energy E, we are sure that the arcs with the same energy E, but different classes, 
cannot have the same value of A\ . Then the velocities £g , <p' cannot coincide (see 
system ©). However, we can't state that collision arcs of different classes determine 
a different set of directions at the point C. What we can certainly assure is that if we 
fix an arc of class q £ I, then for each class q' £ I, not necessarily different from q, we 
can always choose a pair of arcs of class q' , which start at C with directions transverse 
to the velocity with which the arc of class q has arrived at C . 

Thus, for any fixed sequence {(ftj-fcez, qu € /, we can construct infinite collision 
chains {"fk}kez, such that any 7^ is a piece of a periodic orbit of the 2-centre problem 
(L ) of class qu- 

The assumptions of Theorem [T] are all satisfied and we can apply this theorem and 
then Theorem [5] to get our main result Theorem [31 which is now definitely proved. □ 
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A Positivity of entropy 



In this appendix we give a proof of the positivity of the topological entropy of the Poincare map 
for the case at hand. We do it by computing the exponential growth rate of the number of periodic 
orbits. 

By local uniqueness, if ("fk )iez ls a periodic collision chain, then the shadowing orbit is also 
periodic. We compute for each positive integer n the number P n of periodic collision chains with 
period n and take the limit 

p = hm sup . 

n— too Tl 

Let q be fixed and consider only the collision arcs of class q. They are exactly four, corresponding 
to the two possible transverse directions at the centre C. We distinguish two cases: 

i) there are exactly two transverse periodic orbits of the unperturbed problem through C (up to 
reverse the direction of motion), and the arcs are given by the entire periodic orbits; 

ii) there is only one periodic orbit of the unperturbed problem through C (up to reverse the 
direction of motion), with a transverse autointersection at C , and the arcs correspond to parts 
of the periodic orbit between two successive passages through C. 

In case i), for a sequence (7fe i )igz to be a collision chain we must have that 7%. (0) = 7%. (Tj-) ^ 
±7fc 4 (0). Then we cannot have periodic orbits with odd period n = 2m+l, m 6 Z+ and P2 m +l = 
for each m £ Z + . The number of periodic orbits with even period is P^m = 2 2m+1 instead, then 
p = log 2. 

In case ii), for (7fc 4 )iez to be a collision chain, we must have that 7^(0) = ±7fc i+1 (0), because 
there are only four possible velocities at C, divided in two pairs each containing the parallel ones, and 
any arc 7*. arrives at C with velocity 7fc i (Tj. i ) transverse to the initial one 7fe i (0). Then P n = 2 n+1 , 
for any n 6 Z+ , and again p = log 2. 

We conclude that the entropy corresponding to any choice of / C Q + is /it op > log 2. 



References 

[1] S. Aubry, R.S. MacKay, C. Baesens, Equivalence of uniform hyperbolicity for symplectic twist 
maps and phonon gap for Frenkel- Kontorova models, Physica D 56 (1992), 123-134. 

[2] S.V. Bolotin, Nonintegrability of the n-center problem for n > 2, Mosc. Univ. Mech. Bull. 39, 
No. 3, 24-28 (1984); translation from Vestnik Mosk. Gos. Univ., Ser. I, math. mekh. 3 (1984), 
65-68. 

[3] S.V. Bolotin, Influence of singularities of the potential energy on the integrability of dynamical 
systems, J. of Applied Math, and Mech. 48 (1985) No. 3, 255-260. 

[4] S.V. Bolotin and R.S. MacKay, Periodic and Chaotic Trajectories of the Second Species for the 
n-Centre Problem, Celest. Mech. & Dyn. Astr. 77 (2000), 49-75. 

[5] S.V. Bolotin and R.S. MacKay, Non-planar second species periodic and chaotic trajectories for 
the circular restricted three-body problem, Celest. Mech. & Dyn. Astr. 94 (2006), No. 4, 433-449. 

[6] S.V. Bolotin and P. Negrini, Regularization and topological entropy for the spatial n-center 
problem, Ergod. Th. & Dynam. Sys. 21 (2001), 383-399. 

[7] S.V. Bolotin and P. Negrini, Chaotic behavior in the 3-center problem, J. Differential Equations 
190 (2003), 539-558. 

[8] C.L. Charlier, Die Mechanik des Himmels, Bd. I, II, Vcrlag Von Veit & Comp., Leipzig, 1902. 

[9] Y. Duan and J. Yuan, Periodic orbits of the hydrogen molecular ion, Eur. Phys. J. D 6, (1999), 
319-326. 

[10] Y. Duan, J. Yuan and C. Bao, Periodic orbits of the hydrogen molecular ion and their quanti- 
zation, Phys. Rev. A 52, 5 (1995), 3497-3502. 

[11] J. Font, A. Nuncs, C. Simo, Consecutive quasi collisions in the planar circular RTBP, Nonlin- 
earity 15 (2002), 115-142. 

[12] J. Font, A. Nuncs, C. Simo, A numerical study of the orbits of second species of the planar 
circular RTBP, Celest. Mech. & Dyn. Astr. 103 (2009), 143-162. 

[13] A. Knauf and I. A. Taimanov, On the integrability of the n-centre problem, Math. Ann. 331 
(2005), 631-649. 



22 



